



THE GEOMETRY OF INTEGRABLE AND 
SUPERINTEGRABLE SYSTEMS 
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Abstract. The group of automorphisms of the geometry of an 
integrable system is considered. The geometrical structure used 
to obtain it is provided by a normal form representation of intc- 
£*\j . grable systems that do not depend on any additional geometrical 

» . | structure like symplectic, Poisson, etc. Such geometrical structure 

provides a generalized toroidal bundle on the carrier space of the 
system. Non-canonical diffcomorphisms of such structure generate 
alternative Hamiltonian structures for complete integrable Hamil- 
tonian systems. The energy-period theorem provides the first non- 
trivial obstruction for the equivalence of integrable systems. 
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1. Introduction 



After the Erlangen Programme written by Felix Klein with the col- 
laboration of Sophus Lie (Lie visited Klein for two months, just before 
the programme was written) it is by now clear that a "geometry" or 
a "geometrical structure" on a manifold M amounts to select a sub- 
group of the group of diffeomorphisms of M |K172] . |K192] (a recent 

VO \ transcription is available on the arXiv: 0807.3161 vl [math. HO]). 

In Physics we are quite familiar with this correspondence, for in- 
stance we have the Poincare group for the Minkowski space-time in 
Special Relativity, the group Diff(IR 4 ) in general relativity, the group 
of symplectomorphisms in Hamiltonian dynamics, contact transforma- 
tions, the unitary group in Quantum Mechanics and so on. In all 
previous cases, the subgroup of the group of diffeomorphisms deter- 

^ mined by the geometry characterizes the geometry itself, for instance, 

if (Mi,Ui) and (M 2 ,uj2) are two symplectic manifolds and there is a 
group isomorphism $: Sp(Mi,u>i) — > Sp(M 2 ,o;2) of the corresponding 
groups of canonical transformations, i.e., symplectic diffeomorphisms 
defined by each one, then the two symplectic manifolds are (up to 
a conformal constant) symplectically equivalent, i.e., there exists a 
diffeomorphisms cp: Mi — y M 2 such that (p*u 2 = cuj\ [Ba86j-[Ba88j. 
Similar results were established by Grabowski for Poisson and Jacobi 
manifolds |Gr00j . Thus the autormorphisms determined by some geo- 
metrical structures are essentially inner as it happens for the group of 
unitary transformations of a Hilbert space. 



^n leave of absence from the Depto. di Fisica Teorica, Univ. Federico II di 
Napoli, Via Cintia, Napoli, Italia. 
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In this paper we would like to identify the appropriate subgroup of 
Diff(M) associated with integrable and superintegrable systems. 

Thus if we have an integrable Hamiltonian dynamical system T de- 
fined on a symplectic manifold (M, u) the naive thought that the sub- 
group of the group of diffeomorphisms determined by it should be a 
subgroup of the group of canonical transformations of u> is immedi- 
ately shown to be inadequate because integrable ssytems always admit 
alternative Hamiltonian descriptions and we would not know which 
canonical transformations to consider. 

To show how a geometrical structure is characterized by a tenso- 
rial object (and the associated invariance subgroup) we consider a few 
examples. We will consider first a linear structure on a manifold. A 
possible linear structure on a manifold M is characterized by a com- 
plete vector field A such that: 

(1) There exists only one point, say mo G M, such that A (mo) = 0. 

(2) The eigenvalue problem: £&f = • /, / G J-{M) has only 
trivial solutions, i..e, / = constant if M is connected. 

(3) The eigenvalue problem: £&f — f, f G J-'(M) has dimM = n 
independent solutions fk such that dfi A • • • A df n ^ 0. 

Let us remark that the completeness condition on A allows to "gene- 
rate" all of M starting with a properly chosen transversal codimen- 
sion one submanifold in a neighborhood of itlq considered as a set of 
"initial conditions". Any such vector field A identifies a subgroup of 
Diff(M) by requiring <p*(A) = A. Thus the subgroup GL(M,A) of 
Diff(M) of diffeomorphisms tp preserving A is exactly the group of lin- 
ear isomorphisms GL(n, M) with n = dimM where we use the global 
chart provided by the functions fj to identify M with M. n . In infin- 
itesimal terms, linear vector fields are then solutions of [X, A] = 0. 
Moreover if we have two linear structures (Mi, Ai) and (M 2 , A 2 ) de- 
fined in two manifolds Mi and M 2 , and there is a group isomorphisms 
\l/: GL(Mi,Ai) — > GL(M 2 ,A 2 ), then there exists a diffeomorphism 
ip: Mi -^ M 2 such that V*Ai = A 2 . 

The notion of linear structure can be weakened by replacing con- 
ditions (l)-(3) above by asking that the set of points in M satisfying 
A(ra) = define a smooth submanifold Z of M of dimension k. Condi- 
tion (2) will be rephrased by asking that functions / such that £&f = 0, 
defines an Abelian algebra whose spectrum is diffeomorphic to Z and, 
finally, condition (3) will be substituted by demanding that C&f = f, 
define n — k functionally independent fibrewise-linear functions. Such 
"partial" linear structure is actually equivalent to a vector bundle struc- 
ture on the manifold M and the corresponding vector bundle autor- 
mophisms are selected by the condition Lp*A = A. 

Another example in the same vein is concerned with possible tangent 
and cotangent bundle structures on a manifold M. 
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We recall that a tangent bundle structure on M is identified by a 
pair (A,S), where A defines a partial linear structure on M, i.e., a 
vector bundle structure, and S is a (1, l)-tensor field (see for instance 
[Figgj ) such that kerS = ImS (that implies S 2 = 0), S(A) = and 
dg — where the "twisted" differential ds is defined by 

(1) dsf(X) = df(S(x)). 

The vector field A is required to define a vector bundle structure on 
M, i.e., to satisfy the modified conditions (1), (2) and (3) above with 
2k = dimM. The vector bundle structure identified by functions in (3) 
becomes the tangent bundle structure of a manifold Q, so M = TQ. 
Notice that the (1, l)-tensor 5* is the soldering form of TQ. In natural 
bundle coordinates (q,v), these tensor fields take the form: 

„ . d . d 

S = dq®—, A = v—. 

ov ov 

The subgroup of Diff(M) identified by ip*A = A and ip*S = S is 
the group of tangent bundle automorphisms. A single manifold can 
be equipped with alternative tangent bundle structures by considering 
diffeomorphisms (p: M — >• M such that tp*S = S but not preserving 
A, we consider A x and A 2 given as Ai = A and A 2 = y?*A. Then M 
becomes a double vector bundle if [A 1; A 2 ] = and M will carry two 
tangent bundle structures TQ\ and TQ 2 . 

Similarly one may define a cotangent bundle structure by means of 
a pair (A, 9) where again, A is a partial linear structure on M, 9 is 
a particular one-form such that d9 is a symplectic form, i^d9 = 9, 
and requiring that solutions of £&f = are pairwise in involution and 
define a maximal Abelian subalgebra with respect to the the Poisson 
bracket defined by u = d9. Notice that this is equivalent to asking 
that the submanifold Q defined by A(m) = is Lagrangian with re- 
spect to d9. Thus M = T*Q and the canonical Liouville one-form of 
T*Q is just the one-form 9 above. Again as in the case of tangent 
bundle structures, alternative cotangent bundle structures on M can 
be constructed by choosing diffeomorphisms if such that cp*9 ^ 9. 

Thus we will try to study the geometry of an integrable system T by 
determining its associated subgroup of the group of diffeomorphisms 
of the manifold by writing the system in "normal form" instead of de- 
terming such a subgroup from the subgroup of diffeomorphisms deter- 
mined by some geometrical structure determined by it. In particular we 
shall consider systems with orbits possessing a compact closure, even 
though systems possessing unbounded orbits are relevant in scattering 
theory for instance. 

To indentify in such a manner the subgroups of diffeomorphisms cor- 
responding to (super) integrable systems let us first recall the standard 
definition for integrable systems. 
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2. The normal form of an integrable system 

2.1. Integrability and alternative Hamiltonian descriptions. On 

a 2r7,-dimensional symplectic manifold (M, u) a vector field T such that 
irOJ = dH is said to be integrable if there exist n functionally indepen- 
dent first integrals fi, . . . , f n , rf/i A • ■ • A df n ^ 0, such that: 

(2) {/i,/ fe } = 0, {#,/ fc } = 0, Vj,£; = l,2,...,n. 

The independence condition dfx A- • ■ A df n ^ may just hold on some 
open dense submanifold of M, or in a weaker form, on some open 
invariant submanifold. We will find this situation for instance when 
considering scattering problems were we prefer to remove closed regions 
composed of bounded orbits from phase space. 

When the system T possesses more than n first integrals, the system 
is said to be superintegrable. In this case, we have /i, . . . , f n +k, dfi A 
■ ■ ■ A df n+ k 7^ 0, and {H, fj} = 0, for all j = 1, 2, . . . , n + k. In the 
particular case of n + k = 2n — 1 the system is said to be maximally 
superintegrable. We remark that most of the known integrable systems 
are also maximally superintegrable. 

To study the subgroup of diffeomorphisms of M appropriate for an 
integrable system, it is convenient to have a "normal form". 

As it was discussed in the introduction, in searching for normal 
forms, it is quite natural to ask which transformations are allowed 
on the sytem. Thus for Hamiltonian systems, it would be natural to 
use canonical transformations, this is to consider transformations of 
the systems belonging to the closed subgroup of symplectic diffeomor- 
phisms of the group of diffeomorphisms of our manifold M. However, 
any integrable system admits alternative Hamiltonian descriptions, i.e., 
there are (u a , H a ), a = 1,2, such that 

ircoi = dHi, 2^2 = dH 2 . 

In this case, which canonical transformations should we use, symplectic 
diffeomorphisms with respect to U\ or with respect to uj<p. 

In this respect it is amusing to recall a quarrel between Levi-Civita 
and Birkhoff around this. Indeed in his 1934 paper "A general survey of 
the theory of adiabatic invariants" [Le34] Levi-Civita feels the need to 
write a section entitled "Birkhoff 's severity against canonical variables 
and methods. Apology for a milder attitude" (p. 430). From the 
point of view of Birkhoff, one should consider the orbit of vector fields 
obtained by acting on T with Diff(M). All the vector fields in the same 
orbit share the same properties. Thus, it makes sense to restrict the 
attention to just one representative. Let us elaborate on this point in 
the particular case of linear systems. 

We notice that a Hamiltonian vector field admits a factorization in 
terms of a Poisson tensor and an exact one-form: 

(3) T = A(dH). 
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Clearly, any vector field in the orbit of T will also be decomposable 
as above. Moreover, if we consider diffeomorphisms if G Diff (M) such 
that ip*r = T, and apply it to the decomposition above ([3]), we find: 

We conclude that any non-canonical transformation which is a sym- 
metry for T, produces alternative Hamiltonian descriptions. In infin- 
itesimal form: if X is an infinitesimal symmetry for the dynamics T, 
[T,X] = 0, we get: 

Cx{iru) = dC x H, 
then 

i r {C x u) = d(C x H), 

and the previous equation provides an alternative Hamiltonian descrip- 
tion for T if £ju is nondegenerate. In relation with alternative descrip- 
tions, we notice that there are additional ways to generate alternative 
Hamiltonian descriptions. Let us consider for instance a (1, l)-tensor 
T such that: 

(4) C T T = 0. 

Let us define as before eq. (pQ), the twisted differential: 

d T f(X) = df(TX). 

Then from any constant of the motion F we will have the closed 2-form: 

ujt,f = ddxF, 

now if ut,f is nondegenerate, it provides an alternative Hamiltonian 
description for T. In fact notice that: 

itojt,f = C-vdrF — d(irdrF), 

but CrdrF = because of the invariance condition (J3J) and the fact 
that F is a constant of the motion. Hence —dF(TT) is a Hamiltonian 
function for T with respect to the symplectic structure ujt,f- 

2.2. Integrability and normal forms. Because of the previous dis- 
cussion we should expect, therefore, that for vector fields with a large 
group of symmetries we will find always alternative Hamiltonian de- 
scriptions. However we should stress that there are alternative Hamil- 
tonian descriptions which are not generated by diffeomorphisms. 

In conclusion, we should accept any diffeomorphism to reduce a vec- 
tor field T to its normal form. Thus, if the orbit through T contains a 
completely integrable system we can concentrate our attention on the 
standard normal form we are familiar with when we construct action- 
angle variables, i.e., we could consider the form: 

r = ^(/)/-. 

dcpj 
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We should therefore study a normal form for integrable systems as 
emerging from the following conditions: 

Definition 1. Given a vector field V we will say that F — ^(fi, ■ ■ ■ , fn)Xj 
is a normal form for T (and that T is integrable) if 

i. There exist n functionally independent first integrals f\,..., f n , 
such that dfi A • • • df n 7^ . 

ii. There exist n complete vector fields Xi, . . . ,X n pairwise com- 
muting [Xj,Xi\ = and independent Xi A • • • X n ^ 0, and 

hi. C Xj fi = for all j,l = l,...,n. 

We should notice that we have dropped the requirement for T to be 
Hamiltonian and consequently that the "frequencies" v^ are derivatives 
of the Hamiltonian function. 

The usual Liouville-Arnold's theorem becomes now a particular way 
to find functions (coordinates) which reduce T to normal form. 

A few remarks are in order here. 

(1) All integrable systems have the same normal form, then what 
distinguishes one system from another if any such distinction 
exists? 

(2) Which aspects of the normal form above for a given integrable 
system are able to discriminate integrable from superintegrable 
systems? 

In connection with the first query, we immediately notice that many 
interesting aspects on the qualitative structure of the orbits of the sys- 
tem are to be extracted from the normal form because we know that 
usually specific integrable systems need not be diffeomorphic among 
them. We may rephrase our questions by investigating how many dif- 
ferent orbits exist in X(M) under the diffeomorphism group when each 
orbit is required to contain at least one element which is completely 
integrable. 

To have an idea of the variety of situations we might be facing we 
shall investigate some variations on the theme of harmonic oscillators. 

2.3. Hamiltonian linear systems. We may consider as a particular 
instance of the analysis performed in the previous sections the isotropic 
harmonic oscillator. Thus, let us consider M = M. 2n = C n and the linear 
system T defined in cartesian coordinates (xj,yj) by: 



d d 

then, 



5 x ' = k, »' «» = -^ 



r = I> (»'£-*<£) 
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Introducing complex coordinates Zj = yj + ixj, we obtain: 
d d _ 

— Zj = ILJZn, —Zj = —lOJZj, 

dt J J ' dt 3 3 

The algebra of first integrals is generated by the quadratic forms Z\Z m 
and because T is proportional to the linear vector field defined by the 
complex structure on C n , we conclude that its group of linear symme- 
tries is GL(n,C). For a given factorization T = A(dH), the homoge- 
neous space GL(n,C)/GL(n,C) fl Sp(2n, R) parametrizes alternative 
Hamiltonian descriptions, however not all of them. 

What will happen then for a generic linear sytems? Given a generic 
linear system, represented by the matrix A, it has a decomposition: 

(5) A = A-H, 

with A a nondegenerate skew-symmetric matrix and H a symmetric 
matrix if and only if Tr/l 2fc+1 =0, k = 0,1,2, .. .. The proof is easily 
given in one direction. Indeed, A T = —HA, and therefore all odd 
powers must have the same property, implying the assertion. 

For nongeneric matrices conditions guaranteeing the existence of the 
factorization ([5]) are more cumbersome and we refer to the paper [Gi98j 
for a full discussion. 

If T is a linear transformation such that TAT -1 = A, we have: 

A = (TAT^HT-yHT- 1 ), 

providing an alternative decomposition of A if T is not canonical. 

For instance, when A is generic all symmetries are generated by pow- 
ers of A. The non-canonical ones are given by even powers, therefore: 

T = e XA2 \ k = l,2,..., 

will be a non-canonical symmetry for any value of A. 

For instance, two alternative descriptions for the isotropic harmonic 
oscillator in M 4 are given by: 

dpi dq x dp 2 oq 2 2 

a d d d d TT . . 

A 2 = tt— A- V —— A — , H 2 = u{p 1 p 2 + qiq 2 ). 

op i oq 2 op 2 dqi 

A particular invariant (1, l)-tensor field is defined by: 

m 9 d , d , d 

T = dqi® h dq 2 <g> h dp x <g> V dp 2 <g> — - . 

oq 2 dqi ap 2 apx 

We may then consider the 2-form ddxF with F — \{p\+p 2 + qi+ q 2 ) 2 
and get: 

dd T F = d(p\ +p\ + q\ + q 2 ) A d(p 1 p 2 + q x q 2 ) + 

+2{pl +p 2 2 + q\ + q 2 2 )(dq 2 A dq x + d Pl A dp 2 ). 
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We finally remark that the selection of a specific decomposition with 
the Hamiltonian being positive definite gives a group of canonical sym- 
metries which is the unitary group: 

GL(n,C) (1 0(2n,R) = U(n), 

therefore the system may be thought of as a "quantum-like" system 
[ErTU] . 

3. The group of diffeomorphisms of an integrable system 

If T is an integrable system possessing a normal form like in Def. [TJ 
then because X 1; . . . , X n are complete and pairwise commuting we can 
define an action of the Abelian group M. n onto M. Moreover, because 
£x J z /fc = 0, we could redefine the vector field Xj to Yj = v^Xj (no 
summation on j) and we would still have pairwise commuting vector 
fields [Yj, Yfc] = for all j, k — 1, . . . , n. Notice that the completeness 
condition will not be spoiled by redefining the vector fields in this 
form, i.e., the vector fields Y& will be complete (but with a different 
parametrization) and we would define an alternative action of W 1 on 
M. 

Thus what seems to matter is not the particular action of M. n but 
rather the integral leaves of the involutive distribution generated by 
Xi, . . . ,X n . In those cases where the leaves are compact, say tori, we 
could require the choice of an action of M. n that factors to an action of 
T n = R n /Z n . Moreover we could select a particular basis of vector fields 
such that Xi, . . . ,X n each generate the action of a closed subgroup. 
With these particular prescriptions we decompose our vector field I\ 
Let us denote the particularly chosen generators of closed subgroups 
by Z 1 ,...,Z n , then: 

T = co j Z r 

Thus when the closure of generic orbits of T are n-dimensional, the 
system does not have additional first integrals. If for some leaves the 
closure of the orbits does not coincide with the full torus, there are 
additional invariant relations. When for any initial conditon, the clo- 
sure of the orbits is some k < n dimensional torus then the system has 
additional first integrals and it is superintegrable. 

When the closure is one dimensional for all initial conditions, the 
system is maximally superintegrable. 

It is now clear that the geometry of integrable and superintegrable 
systems is associated with a toroidal generalized bundle, i.e., projec- 
tions on the manifold M which have fibers diffeomorphic to tori of 
dimension going from one to n. This is exactly the situation that hap- 
pens if M is a compact 2n-dimensional symplectic manifold and the 
dynamics is invariant under the action of a torus group T". Then be- 
cause of Atiyah's convexity theorem, the momentum map J : M — > i* 
is a surjective map onto a convex polytope Pet*. The fibers of J 
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are invariant tori that on points on the interior of P are n-dimensional. 
The fibers corresponding to the boundary of the polytope are lower 
dimensional tori [At82j. 

Then the associated subgroup of diffeomorphisms of M determined 
by T is the subgroup of bundle automorphisms of such generalized 
toroidal bundle. As we will see later on, in connection with specific 
integrable or superintegrable systems, the most important obstruction 
to their being diffeomorphic is the energy-period theorem, that actually 
puts a restriction on the nature of the toroidal bundle of the system 
and in consequence on its group of diffeomorphisms. 

4. Oscillators and nonlinear oscillators 

To illustrate the previous situation we give now a few examples. 

We may consider the isotropic Harmonic oscillator with two degrees 
of freedom. Say M = R 4 , to = J2 a dPa A dq a , and H = \J2 a (Pa + 
ql). In this case on R 4 , = R 4 — {0}, an open dense submanifold, the 
dynamics generates orbits with coincide with their closure, they are 
one-dimensional. The toroidal bundle is provided by: 

S 1 ->• R 4 ->• S 2 x R+. 

We have Rq = S 3 x R + and the dynamics induces the Hopf fibra- 
tion S 1 — y S 3 — y S 2 . The subgroup of diffeomorphisms of R 4 , which is 
selected by the fibration is the group of projectable diffeomorphisms. 
Clearly this large group of symmetries, when applied to a chosen Hamil- 
tonian description will generate many more alternative descriptions. 

It should be remarked that the alternative Hamiltonian description 
provided by dp\ A dq\ — dp 2 A dq 2 , and 

H = \(p\ + ql)-\(pl + ql) 

cannot be derived from the standard one with the positive definite 
Hamiltonian H Q because the diffeomorphism would preserve the signa- 
ture of H (because it cannot map compact energy levels, the ones of 
the Hamiltonian H , into non-compact ones, the ones of H). 

The system is actually superintegrable, indeed the quotient manifold 
under the action of the dynamics is three dimensional instead of two 
dimensional. 

This example generalizes to any finite dimension and we have again 
for the n-dimensional isotropic harmonic oscillator the fibration: 

S l ->■ R^ n -»■ CP"- 1 x R + . 

Again the symmetry group for the dynamics is the group of dif- 
feomorphisms projectable under the previous projection, hence it is 
diffeomorphic to the central extension of Diff(CP n_1 x R + ) by U(l). 
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The one-form ^ J^ fc ( Pk JfTf k i\ k ) = dr has the property that irdr = 

1. Any closed two-form on the quotient manifold and any function on 
the quotient manifold without critical points on the invariant open 
dense submanifold specified by p\ + q\ ^ 0, k — 1, 2, . . . , n, give rise to 
an alternative Hamiltonian desription. 

This example is a normal form for maximally superintegrable systems 
with one-dimensional closed orbits and constant period. 

In higher dimensions we can consider the Hamiltonian H = ^ a u a H a , 
with H a = |(Pq+9„), w a G IKL The subgroup associated with the Hamil- 
tonian vector fields r a , ir a co = dH a , are closed subgroups. When all 
frequencies are pairwise irrational, i.e., oo a /uJb is irrational, the closure 
of the generic orbit of T is the full torus, in this case there cannot be 
additional constants of the motion. When some of the frequencies are 
pairwise rational, the closure of a generic orbit is a torus of lower di- 
mensions. A particular example where the closure of the orbits goes 
from a one-dimensional torus to an n-dimensional one, depending on 
the initial conditioins, is provided by: 



\2 



a 

In this case we may also find invariant relations for particular values 
of the initial conditions. 

This example gives rise to the so called nonlinear oscillators and has 
been considered in quantum mechanics to give interesting consequences 
at the level of Planck's distribution law and for alternative commuta- 
tion relations |Ma97j . [Lo97j . 

5. Obstructions to the equivalence of integrable systems 

If we inquire about the obstructions to the existence of diffeomor- 
phisms conjugating two integrable systems, the energy-period theorem 
provides the first and most important one. We refer to the literature 
[Ne02] , |Go69] for various versions of this theorem, however a "cheap" 
proof may be given as follows. On the carrier space of our dynamical 
system let us consider the two-form dpk A dq k — dH A dt. Consider 
now the evolution t *-¥ t + r, Pk(t) = Pk(t + r ) 5 1 k if) — q k (t + t) and 
H(t) = H(t') for all t' G IKL Because the evolution is canonical, we find 
dH A dr = on the submanifold on which the period is a differentiable 
function on (q h ,Pk)- It follows that the period r and the Hamiltonian 
H are functionally dependent. 

It is well-known that the period of a dynamical system is an at- 
tribute of the vector field which is invariant under diffeomorphisms. It 
follows that if two Hamiltonian systems have different periods, they 
cannot be connected via diffeomorphisms. For instance the isotropic 
Harmonic oscillator and the Kepler problem cannot be connected by a 
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diffeomorphism. Indeed the map connecting solutions of the Harmonic 
oscillator and those of the Kepler problem, the Kustaanhneimo-Stiefel 
map, is a map defined on each energy level, i.e., for those orbits that 
have all the same period. The map changes from one energy-level to 
another (see for instance [Av05] . [Av05bj ). 

A more simple example is provided by H = |(j> 2 + q 2 ) and H' = 
(p 2 + q 2 ) 2 - For the first one, the frequency is independent from the 
energy, while for the second one it depends on the initial conditions. 
The two systems cannot be diffeomorphic. This circumstance was the 
main motivation to introduce the classification of dynamical systems 
up to parametrization, i.e., up to conformal factors [Ib94j. 

We hope we have made clear that the "geometrical picture" we have 
derived is the best one can do because each individual integrable system 
will give rise to infinitely many different situations which cannot be 
classified in a meaningful way otherwise (i.e., identifying a finite or a 
countable family of equivalence classes). 
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